Effective attraction induced by repulsive interaction in a spin-transfer system 
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In magnetic systems with dominating easy-plane anisotropy the magnetization can be described 
by an effective one dimensional equation for the in-plane angle. Re-deriving this equation in the 
presence of spin-transfer torques, we obtain a description that allows for a more intuitive under- 
standing of spintronic devices' operation and can serve as a tool for finding new dynamic regimes. 
A surprising prediction is obtained for a planar "spin-flip transistor" : an unstable equilibrium point 
can be stabilized by a current induced torque that further repels the system from that point. Stabi- 
lization by repulsion happens due to the presence of dissipative environment and requires a Gilbert 
damping constant that is large enough to ensure overdamped dynamics at zero current. 

PACS numbers: 72.25.Pn, 72.25.Mk, 85.75.-d 



In physics, there are cases where due to the presence 
of complex environment a repulsive force can lead to ac- 
tual attraction of the entities. A well known example 
is a superconductor, where the Cooper pairs are formed 
from electrons repelled by the Coulomb forces due to the 
dynamical elastic environment. Here we report a phe- 
nomena of effective attraction induced by the repulsive 
spin-transfer torque in the presence of highly dissipative 
environment. The spin-transfer effect producing the re- 
pulsive torque is a non-equilibrium interaction that arises 
when a current of electrons flows through a non-collinear 
magnetic texture [H, 0, 01 • This interaction can become 
significant in nanoscopic magnets and is nowadays stud- 
ied experimentally in a variety of systems. Its manifesta- 
tions - either current induced magnetic switching [J] or 
magnetic domain wall motion [5fl - serve as an underlying 
mechanism for a number of suggested memory and logic 
applications. 

Here we consider a conventional spin-transfer device 
consisting of a a magnetic polarizer (fixed layer) and a 
small magnet (free layer) with electric current flowing 
from one to another (Fig 1). Both layers can be de- 
scribed by a macro-spin model due to large exchange 
stiffness. The free layer is influenced by the spin transfer 
torque, while the polarizer is too large to feel it. Magnetic 
dynamics of the free layer is described by the Landau- 
Lifshitz- Gilbert (LLC) equation with the spin transfer- 
torque term 0, Q . 

The solutions of LLG are easy to find for the simplest 
easy axis magnetic anisotropy of the free layer. There 
exists a critical current at which the free layer either 
switches between the two minima of magnetic energy, 
or goes into a state of permanent precession, powered by 
the current source [2, H, 0]. The same basic processes 
happen in the case of realistic anisotropics, however the 
complexity of the calculations increases substantially. In 
a nanopillar device Q one additionally finds that stabi- 
lization of magnetic energy maxima is possible ( "canted 



states" [fij) and that multiple precession modes exist with 
transitions between them happening as the current is in- 
creased 0, !, E3- The anisotropy of a nanopillar device 
is a combination of a magnetic easy plane and magnetic 
easy axis directed in that plane. Experimentally, the 
easy plane anisotropy energy is usually much larger than 
the easy axis energy, i.e. the system is in the regime 
of a planar spintronic device [IH (Fig. [T]) . This limit of 
dominating easy plane energy is characterized by another 
simplification of the dynamic equations fl2l Il3j |. which 
comes not from the high symmetry of the problem, but 
from the existence of a small parameter: the ratio of the 
energy modulation within the plane to the easy plane en- 
ergy. The deviation of the magnetization from the plane 
becomes small, making the motion effectively one dimen- 
sional. 

In this paper we present a general form of effective 
planar equation describing a macrospin free layer in the 
presence of spin transfer torques. Its relationship to the 
first order expansion in the current magnitude used in 
Ref. [HI is discussed at the end. We then use this equa- 
tion to study the "spin-flip transistor" : a planar device in 
which the spin polarizer is perpendicular to the direction 
favored by the magnetic anisotropy energy. It was pre- 
dicted [14J that the competition between the anisotropy 
and spin transfer torques leads to a 90 degrees jump of 




FIG. 1: Planar spin-transfer devices. Hashed parts of the 
devices are ferromagnetic, white parts are made from a non- 
magnetic metal. 
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the magnetization at the critical current. Whether the 
jump happens into the parallel or antiparallel state with 
respect to the polarizer is determined by the direction of 
the current. 

Here it is shown that the behavior of the spin- flip tran- 
sistor is more complicated than expected from the simple 
picture above. Namely, the current inducing a jump into 
the parallel direction can also stabilize the antiparallel 
direction. This conclusion is certainly counter-intuitive 
because the spin torque repels the magnetization from 
this already unstable saddle point of the energy. How- 
ever, a combination of two destabilizing torques manages 
to result in a stable equilibrium. We will see that this 
happens due to the dissipation terms and a sufficiently 
large (but still small compared to unity) Gilbert damping 
constant is required to observe the phenomena. 

The magnetization of the free layer M = Mn has a 
constant absolute value M and a direction given by a 
unit vector n(<). The LLG equation 0, 0| reads: 
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Here 7 is the gyromagnetic ratio, E(n) is the magnetic 
energy of the free layer, and a is the Gilbert damping 
constant. The second term on the right is the spin trans- 
fer torque, where s is a unit vector along the direction 
of the polarizer, and the spin transfer strength it(n) is 
proportional to the electric current / @, [H} . In general, 
spin transfer strength is a function of the angle between 
the polarizer and the free layer u(n) — /[(n • s)] I, with 
the function /[(n • s)] being material and device specific. 
Equation {T]) can be written in polar angles (9(t), cf>(t)): 
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with tangent vectors = [z x n]/sin6>, eg = [e^ x n]. 
The easy plane is chosen at 9 — tt/2, and the mag- 
netic energy has the form E = (K±/2) cos 2 9 + E r {9, (ft), 
where E r is the "residual" energy. In the planar limit, 
K± — ► 00, the energy minima are very close to the easy 
plane and the low energy solutions of LLG have the prop- 
erty 9{t) = vr/2 + 59 with 59 0. Equations © can 
then be expanded in small parameters \E r \/K± <C 1, 
^(n^/i^x <C 1. Assuming time-independent u and s we 
obtain an effective equation of the in-plane motion 

1 v ; 7 dE e rr 

which has has the form of the Newton's equation of mo- 
tion for a particle in external potential E e ff((ft) with a 
variable viscous friction coefficient a e // {(ft)- The expres- 
sions for the effective friction and energy are 

a e//(0) = a - (r^ + T e )/u±_ , (4) 
T = (dF^/d(ft) d=v/2 , T e = (dF g /d9) g=7l/2 , 



(3) 



E e f f (<f>)=E r (w/2,<p) + AE(cl>) 
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Equation §5§ with definitions (|4l5p gives a general de- 
scription of a planar device in the presence of spin trans- 
fer torque. At non-zero current the effective friction can 
become negative (see below), and the effective energy is 
not necessarily periodic in (ft (e.g. in the case of "mag- 



netic fan" 13 



15j). Physically this reflects the possibility 
of extracting energy from the current source, and thus 
developing a "negative dissipation" in the system. 

In many planar devices the polarizer direction s lies 
in the easy plane, 9 S = tt/2, with a direction defined 
by the azimuthal angle <p s . At the same time the resid- 



ual energy has a property (8E r /d9) 
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not shift the energy minima away from the plane. We 
will also use the simplest form /[(n • s)] = const for the 
spin transfer strength. A more realistic function will not 
change the result qualitatively and can be easily used if 
needed. With these restrictions the effective friction and 
the energy correction get the form: 
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In a spin-flip transistor the polarizer direction is given 
by <fi s — tt/2. Following Ref. [3, we consider in-plane 



anisotropy energy E r (ir/2,i 



- (if 1 1/2) cos <p corre- 



sponding to an easy axis. Then the effective friction 
is a e ff = a + (2u sin (ft) / lu ± and effective energy equals 
(■y/M)E eff = — - u 2 /uj ± )/2] cos 2 (ft + const with 
W11 = jK\\/M. Equilibrium points (ft = 0, ±7r/2,7r are 
the minima and maxima of the effective energy, and do 
not depend on u. Stability of any equilibrium in one di- 
mension depends on whether it is a minimum or a maxi- 
mum of E e f f and on the sign of a e f / at the equilibrium 
point. It is easy to check, that out of four possibilities 
only an energy minimum with a e f / > is stable. In the 
case of a spin-flip transistor the energy landscape changes 
above a threshold |u| > ^uj\\u>±: the energy minima at 
(ft = 0, 7r become maxima, and, vice versa, the energy 
maxima at (ft = ±7r/2 switch to minima. Effective fric- 
tion at (ft — 0, 7r is positive independent of u, while at 
eft = ±7r/2 it changes sign at u = ^auij_/2. 

The behavior of the spin-flip transistor is summarized 
in a switching diagram Fig. [2] plotted on the plane of the 
material characteristic a and the experimental parame- 
ter u ~ J. For definiteness we will discuss a current with 
u > 0. The effect of the opposite current is completely 
symmetric. For small values of Gilbert damping one ob- 
serves stabilization of the (ft = tt/2 (parallel) equilibrium 
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FIG. 2: Switching diagram of the spin-flip transistor. In each 
zone one or two arrows show the possible stable directions of 
the free layer magnetization. Directions of the easy axis and 
spin polarizer are defined in the right bottom corner. Angular 
dependencies of a e f / and E e f / are given in insets. Stable sub- 
regions "b" and "c" differ in overdamped vs. underdamped 
approach to the equilibrium. 



to which the spin torque attracts the magnetization of 
the free layer, while the opposite (antiparallcl) direction 
remains unstable. This is in accord with the results of 
Ref. [3- However, when the damping constant is larger 
than the critical value a„ = 2yJZJyJZT[i a window of sta- 
bility of the antiparallel equilibrium opens on the dia- 
gram. Since a -C 1, a sufficiently large easy plane energy 
is required to achieve a, < a « 1. 

If one thinks about the stability of the (8, <fi) — 
(tt/2, — 7t/2) equilibrium for u > in terms of Eq. this 
prediction seems completely unexpected. The anisotropy 
torques do not stabilize this equilibrium because it is a 
saddle point of the total magnetic energy E, and the 
added spin transfer torque repels n from this point as 
well. The whole phenomena may be called "stabilization 
by repulsion" . To check the accuracy of the planar ap- 
proximation ([3|), the result was verified using the LLG 
equations ([2|) with no approximations for the axis-and- 
plane energy E = (Kj_/2) cos 2 9 - {K\\/2) sin 2 9 cos 2 <j>. 
Calculating the eigenvalues of the linearized dynamic ma- 
trices § at the equilibrium points (7r/2,±7r/2) we ob- 
tained the same switching diagram and confirmed the 
stabilization of the antiparallel direction. Typical trajec- 
tories n(i) numerically calculated from the LLG equa- 
tion with no approximations are shown in Fig. [3] to illus- 
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FIG. 3: Typical trajectories of n(t) for wm/wj. = 0.01, a = 
f .5 a,. The plot labels correspond to the regions in Fig. [2] the 
current magnitude is given in the units of u/ ' yfu\\Ju)± and we 
look at the stability of the <f> = —tt/2 equilibrium: (a) 0.93, 
unstable (b): 1.08, stabilized with overdamped approach (c): 
1.38, stable, but with oscillatory approach (d): 1.53, unstable; 
a stable cycle is formed around the equilibrium. 



trate the predictions. At it > ^/SJfjoJJ the <f> = — 7r/2 
equilibrium is stabilized. In accord with the predic- 
tions of Eqs. (J!]),©, the wedge of its stability consists 
of two regions (b) and (c) characterized by overdamped 
and underdamped dynamics during the approach to the 
equilibrium. The dividing dashed line is given by u — 
u)\\/ot + It was checked that small deviations 

of the polarizer s from the (tt/2, 7t/2) direction do not 
change the behavior qualitatively. Larger deviations 
eventually destroy the effect, especially the out-of-plane 
deviation which produces the "magnetic fan" effect [l5| 
leading to the full-circle rotation of 4> in the plane. 

As the current is further increased to u > auj±/2, the 
antiparallel state looses stability and the trajectory ap- 
proaches a stable precession cycle (Fig. Eld)). The exis- 
tence of the precession state is easy to understand from 
(|3|) viewed as an equation for a particle in external poten- 
tial. Just above the stability boundary the effective fric- 
tion a e ff(4>) is negative in a small vicinity of <p = —tt/2, 
and positive elsewhere. Within the a e ff < region the 
dissipation is negative and any small deviation from the 
equilibrium initiates growing oscillations. As their am- 
plitude exceeds the size of that region, part of the cycle 
starts to happen with positive dissipation. Eventually 
the amplitude reaches a value at which the energy gain 
during the motion in the a e f / < region is exactly com- 
pensated by the energy loss in the a e f / > region: thus 
a cycle solution emerges. The effective planar description 
allows for the analysis of the further evolution of the cy- 
cle with transitions into different precession modes, which 
will be a subject of another publication. 
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The fact that a > a* condition is required for the 
stabilization means that dissipation terms play a crucial 
role entangling two types of repulsion to produce a net 
attraction to the reversed direction. Note that an in- 
terplay of a strong easy plane anisotropy and dissipa- 
tion terms produces unexpected effects already in con- 
ventional (u = 0) magnetic systems. The effective planar 
equation ([3|) at u = was discussed in Ref. UH It was 
found that the same threshold a* represents a bound- 
ary between the oscillatory and overdamped approaches 
the equilibrium. Above a* the familiar precession of a 
magnetic moment in the anisotropy field is replaced by 
the dissipative motion directed towards the energy mini- 
mum. When the easy plane anisotropy is strong enough 
to ensure a ^> a*, one can drop the second order time 
derivative term in Eq. ([3]) and use the resulting first or- 
der dissipative equation. In the presence of spin transfer, 
a e ff(4>, u) depends on the current and can assume small 
values even for a 3> a* , thus no general statement about 
the 4> term can be made. 

The simplest easy axis energy expression E r (ir/2,cf>) — 

— (K\\/2) cos 2 cj> happens to have the same angular de- 
pendence as AE((j)) given by Eq. (J6]). Due to this spe- 
cial property the energy profile flips upside down at 
u = y/uJjjaJjT. For a generic E r (jr/2, (f) with minima at 
4> = 0, 7r and maxima at (f> = ±7r/2 the nature of equilib- 
ria will change at different current thresholds. This will 
make the switching diagram more complicated, but will 
not affect the stabilization by repulsion phenomena. Sim- 
ilar complications will be introduced by a generic /[(n-s)] 
angular dependence of the spin transfer strength. 

In Ref. [13| the known switching diagram for the 
collinear (<f> s = 0) devices [1, [j| [l(| were reproduced 
by equation (3]) with E eff = E r (-K /2,4>). The AE 
term ([6]) was dropped as being second order in small 
u. This approximation gives a correct result for the 
following reason. In a collinear device (7 /M)E e f / = 

— [(a; 1 1 + u 2 /u!±)/2] cos 2 4> + const and the current never 
changes the nature of the equilibrium from a maximum 
to a minimum. Consequently, dropping AE does not af- 
fect the results. As was already noted in Ref. [3, the first 
order expansion in u is insufficient for the description of 
a spin-flip transistor, where the full form ^ is required. 

In summary, we derived a general form of the effec- 
tive planar equation ([3]) for a macrospin free layer in 
the presence of spin transfer torque produced by a fixed 
spin-polarizer and time- independent current. Qualitative 
understanding of the solutions of planar equation is ob- 
tained by employing the analogy with a one-dimensional 
mechanical motion of a particle with variable friction co- 
efficient in an external potential. The resulting predic- 
tive power is illustrated by the discovery of the stabi- 
lization by repulsion phenomena in the spin-flip device. 
Such stabilization relies on the form of the dissipative 
torques in the LLG equation and happens only for a large 
enough Gilbert damping constant. The new stable state 



and the corresponding precession cycle can be used to 
engineer novel memory or logic devices, and microwave 
nano-generators with tunable frequency. 

To observe the phenomena experimentally, one has to 
fabricate a device with a > a*, and initially set it into a 
parallel or antiparallel state by external magnetic field. 
Then the current is turned on and the field is switched off. 
Both states should be stabilized by a moderate current 
yfu>\\uj± < u < aui±/2, but cannot yet be distinguished 
by their magnetoresistive signals. The difference can be 
observed as the current is increased above the auj±/2 
threshold: the parallel state will remain a stable equilib- 
rium, while the antiparallel state will transform into a 
precession cycle and an oscillating component of magne- 
toresistance will appear. 
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